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ABSTRACT 

The assessment of turbulence model performance in predicting flow fields that are 
directly relevant to industrial needs has become increasingly important. In aerody- 
namics, many problems exist for which the performance of some turbulence models is 
ambiguous at best and incorrect at worst. The application of many different turbu- 
lence models to a particular problem is a lengthy process that can lead to confusing 
results simply as a result of the volume of data produced. Given the motivation 
to develop turbulence models that do have some applicability to solving industry’s 
problems, this study focuses on the prediction of airfoil characteristics, including 
lift and drag, over a range of Reynolds numbers. Two different turbulence models, 
which represent two different types of models, are tested. The first is a standard 
isotropic eddy-viscosity two-equation K — e model, and the second is an explicit al- 
gebraic stress model (EASM). The EASM model is an extension of the K — e model 
because it introduces nonlinear dependency of the Reynolds stresses on the mean 
strain and rotation-rate tensors and, therefore, better accounts for the turbulent 
stress anisotropy more effectively. The turbulent flow field over a general- aviation 
airfoil (GA(W)-2) at three Reynolds numbers ranging from 2.1 x 10 6 to 6.3 x 10 6 is 
studied. Experimental pressure coefficients are compared with model predictions for 
Re = 4.3 x 10 6 . At each Reynolds number, predicted lift and drag values at different 
angles of attack are compared with experimental results, and predicted variations of 
stall locations with Reynolds number are compared with experimental data. Finally, 
the size of the separation zone predicted by each model is analyzed, and correlated 
with the behavior of the lift coefficient near stall. In summary, the EASM model is 
able to predict the lift and drag coefficients over a wider range of angles of attack than 
the K — e model for the three Reynolds numbers studied. However, both models are 
unable to predict the correct lift and drag behavior near the stall angle, and for the 
Re = 2.1 x 10 6 case, the K — £ model did not predict separation on the airfoil near 
stall. 



I INTRODUCTION 


The optimization of the performance and, therefore, the design of wing configu- 
rations is a challenging task that requires efficient use of numerical prediction tech- 
niques. In fact, an accurate predictive technique is a mandatory requisite. The 
purpose of this study is to assess the predictive capabilities of two types of turbulence 
models for the simpler two-dimensional airfoil configurations. The two types of mod- 
els represented are an isotropic eddy-viscosity K — e model and an explicit algebraic 
stress model (EASM). They are applied to the flow over a general-aviation airfoil at 
various angles of attack and Reynolds numbers. 

One problem with this type of study is the lack of detailed information in regard to 
the experimental setup, as well as the lack of necessary information about initial and 
inflow conditions. Nevertheless, some reports exist that provide sufficient information 
so that a study such as the one outlined here can be performed. One such study 1 
details the aerodynamic characteristics of a 13-percent- thick general-aviation airfoil at 
various angles of attack and Reynolds numbers; the results of this study are used here 
to determine the ability of the models to predict some characterizing features of airfoil 
flows, such as the lift and drag coefficients and the pressure coefficient distributions 
at various flow conditions. 

Other studies have examined the performance of K—e and algebraic stress models 
in predicting the aerodynamic characteristics of airfoil flows. Although the exact form 
of the two-equation models and ASM’s differs between these studies (as well as this 
one), the intent is to examine whether the increased physics introduced into the 
ASM have a significant impact on the results. Because the ASM represents stress 
anisotropy effects better, its utilization should give some indication of whether a full 
Reynolds stress closure is worthwhile. Davidson and Rizzi 3 applied a two-equation 
model and an ASM to the ONERA-A airfoil flowfield at a single chord Reynolds 
number Re of 2.1 x 10®. Their results showed that the ASM predicted the variation 
of Ci and the stall location better than the K — e model (and the Baldwin-Lomax 
model, which was also tested). In addition, the calculated pressure and skin-friction 
coefficients agreed very well with the experimental results. 

Recently, Davidson 3 performed a follow-on study with a K — e and a Reynolds 
stress model at the same conditions as the Davidson and Rizzi work. Once again, the 
higher order model outperformed the isotropic eddy-viscosity K — e model. No stall 
was predicted by the K — e model, whereas stall location was accurately predicted 
by the Reynolds stress model, as were the velocity profiles on the airfoil and in the 
wake. 
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Lien and Leschziner 1 * * 4 compared the performance of a K —e, a Reynolds stress, and 
a nonlinear K — e (ASM) model on the ONERA-A airfoil as well, at the same chord 
Reynolds number used in the Davidson studies. Their study included the development 
of an improved ASM and Reynolds stress model to better predict the flow field. Their 
results suggested that the Reynolds stress model outperforms both the K — e and 
nonlinear K — e models in predicting the aerodynamic characteristics of the airfoil. 

In the present study, the comparisons include the effects of Reynolds-number 
variation. In this study, no attempt is made to optimize the models for better pre- 
dictive capabilities; the K — e model of Speziale et al. 5 and an ASM of Gatski and 
Speziale 6 are used. The algebraic stress model 6 is an explicit polynomial represen- 
tation that is discussed in the next section. Comparisons of predicted and mea- 
sured pressure coefficients are shown at two angles of attack at Re = 4.3 x 10 6 . For 
Re = 2.1 x 10 6 , 4.3 x 10 6 , and 6.3 x 10 6 , the variation of lift and drag coefficients Ci 
and Cd is shown as a function of angle of attack a and Ci, respectively. The variation 
of stall angle with Reynolds number is also examined, and the predicted results are 
compared with experimental results. Finally, to better understand the physics behind 
the stall-angle variation predicted by the K — t model and the EASM, comparative 
plots of separation-zone size are shown for the three Reynolds numbers. These results 
provide an indication of the capabilities of the models, and a direction for the next 
step in evaluating model performance and, ultimately, in improving the models. 


II Theoretical Formulation 


The mean quantities of an incompressible turbulent flow are obtained from a 
solution of the Reynolds- averaged Navier-Stokes equations given in nondimensional 
form as 

^ = 0 


dUi _ dtii 

dt Uj dxj 


dij 

dp 


1 doij duty 


(1) 


( 2 ) 


dii Re dij dxj 

where tt,- and p are the mean velocity and pressure, respectively, and Re is the 

Reynolds number based on free-stream conditions and airfoil chord. The right-hand 

side of the equation also contains both the viscous stress tensor <f,j, which is propor- 
tional to the mean strain rate, and the turbulent second-moment correlation tensor 
T{j = ujttj, which requires closure. 

In the incompressible formulation for the mean motion, only the turbulent second- 
moment Reynolds stress tensor needs to be modeled. For the higher order closures, 
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it is best to begin with the transport equation for the T,j correlation: 



dTij , d , n , rrd , d ^ijk , 

dt + dx k {ukTtl) " P,J + n<j e,3+ dx k +l/ dx k dx k 

(3) 

where the right-hand side represents the rate of change of r,j produced by the tur- 
bulent production P,j, the deviatoric part of the pressure strain- rate correlation II fj, 
the turbulent dissipation rate e,j, and the turbulent diffusion D\j k . These terms are 

given by 

= duj dUi 

Pii ~ Tik dx t ~ Tik dx k 

(3a) 



(3b) 


2 c du'i du'j 

c„- 3 et v -2i ’ dxi g Xi 

(3c) 

and 




D' jt = -[«r«K + WAk + «SM 

(3d) 


In the incompressible formulation, the deviatoric part of the dissipation rate is 
generally absorbed into the pressure-strain correlation to account for any anisotropic 
dissipation effects. However, models for the anisotropic dissipation rate can be 
derived that negate this assimilation into the pressure-strain correlation. 7 ’ 8 When the 
deviatoric part of the tensor dissipation rate is assimilated into the pressure-strain 
correlation, the isotropic dissipation rate e is obtained from a modeled transport 
equation. 

The common high-Reynolds-number form of the isotropic solenoidal dissipation- 
rate equation that is used is given by 

d 2 e 


dt + dxj ^ + dx^ tj + "dxjdxj 


with 


and 


V t = -C tl —Tu 


e dUi 
~K Tii ~dxj 


V. = C t2 - 


(4) 


(4a) 


(4b) 


where the first and second terms on the right side are the production and destruction 
of dissipation, respectively, and the last term is the viscous diffusion. The turbulent 
dissipation-rate diffusion term is modeled for the two-equation formulation used in 
this study as 

»< - £ 1 ; < 5 > 
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The closure constants C € i, C € 2 , and a t must be specified to complete the closure. 

This brief outline completes the specification of the transport equations for the 
solution of the second-moment closure model. To this point, no specific models have 
been proposed for terms such as the pressure- strain correlation or turbulent diffusion. 
Many models exist in the literature; the reader is referred to the many review articles 
on the subject. 9-12 

From an engineering standpoint, this second-moment formulation is not without 
risk. Extensive computer resources and sophistication are necessary in implementing 
and numerically solving these equations over a two-equation formulation; furthermore 
no unambiguous assessment has been made of the range of validity of the models, 
which is, of course, unfortunate because these equations represent the inclusion of 
extensive physics and mathematical rigor that should lead to improved predictive 
technology. For this reason, an intermediary level can be introduced into the hierarchy 
of closure models that incorporates the key features of the second-moment closures 
but does not involve more effort than a simple isotropic eddy- viscosity two-equation 
model. 

An effective compromise between the full second-moment closure and the two- 
equation model is the ASM . This approach extracts an algebraic relationship between 
the turbulent Reynolds stress and mean velocity field from the Reynolds stress formu- 
lation just discussed. The models are derivable from equilibrium hypotheses imposed 
on both the convective and diffusive terms. The convective equilibrium hypothesis is 
represented through the stress anisotropy tensor 

, - §*<«) 
b '‘ 2 It 

and is given by 

Dbjj _ Drjj _ Tjj DK _ 0 
Dt ~ Dt K Dt 

where K = r,,/2. The equilibrium hypothesis for the diffusive terms is given by 
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dx k 


D\ ik + v 


\ 

dxkdxk) 


= 0 


( 8 ) 


If we use these constraints, as well as the definition of the stress anisotropy, and 
consider the fact that the mean velocity gradient is the sum of the strain- and rotation- 
rate tensors (Sij and W,- j , respectively), then Equation (7) can be expanded as 


(V — e)bij = — K (bikSjk + bjk^ik — ^bmnSmnSij'j 

- K(bi k Wjk + bjkWik) + \nfj 


(9) 
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where V = Pa/2, ~5ij = ( diii/dxj + dujfdxi)/ 2, and Wij = ( dui/dij — diij/dxi)/ 2. 
This relation is the traditional starting point for the algebraic stress formulation. The 
pressure-strain correlation model Ilfj of Speziale, Sarkar, and Gatski (SSG) 13 is used 
in Equation (9) to close the equation. 

For a two-dimensional flow, the exact integrity basis polynomial representation 6 
of Equation (9) is given by 

tv; = - 2v\ [3* + ^ @ ik W kj + J jk W ki ) 

( 10 ) 

- 2 ^ (5 ik S kj - 


where 


«r = c: 


K 2 


e 

rv 3(1 + if) 

" 3 + 1] 2 + 6CV + 6C 2 


<*i 


v = — (3*3«) 1/a , C = — (W.iWo) 1 / 2 

€ £ 


ai=Q-C 2 )|, a 2 = (2-C 4 )§, a 3 = (2-C 3 )|, g = 
For the SSG model, 13 these coefficients are 


+ C 5 — I 


( 11 ) 

( 12 ) 

(12a) 

(12b) 


c?! = 6.8, Cj = 0.36, C 3 = 1.25, C, = 0.40, C, = 1.88 (13) 


The application of the model to a variety of compressible bump and airfoil flows 14 
has shown that its robustness is diminished by using the C* coefficient given in Equa- 
tion (12). The robustness diminishes because during the initial transient phases and 
in the free-stream region the value of C’ attains values < 0.2ai « 0.023. This low 
value of C’ causes regions of high shear to form to maintain the turbulence, which 
in turn leads to numerical instabilities. To avoid these problems, a modified form of 
Equation (10) is used: 

TV; = \Kf>ij 2v t ~5 j 

K r / 1 \ i (l^) 

- 21/’’ — ja 2 k j — 2a 3 \l>i k !> k j — 

where C* is now given by 15 


3(1 4- rj 2 ) + 0.2(>? 6 + C 6 ) 

M 3 + »7 2 + 6C V + 6C 2 + ri« + C 6 1 


(15) 
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and v ** is now introduced and defined by 


with 




™ 3(1 + I)») 

* 3 + 7j 2 + 6CV + 6C ! + 1 )' + c 6 


<*1 


( 16 ) 


( 17 ) 


For large values of rj and £, C* and C** are limited by 0.2c*i and 0, respectively. 
In near equilibrium conditions, Equation (15) and (17) reduce to the value given in 
Equation (12). 

The nonlinear algebraic constitutive equation given in Equation (14) is then cou- 
pled with an incompressible two-equation formulation. For the K — e formulation, 
the governing equations are subsets of the second-order formulation shown earlier in 
Equation (3) and (4). These governing equations are given by 15 


DK ~ d 

■DF = p-£ + aiT 


(' 


+ 


uu\ dK_ 
(T k J dx k 


and 


with 


= r Lv-c + A 

Dt ° tl K V Cc2 *K + dx k 


„ K 2 

vti = Ctf — 
e 




(18) 

(19) 

(19a) 
(19b) 

(19c) 

« = 0.41 

(C„ - C„) 

where u T is the shear velocity and y is taken as the distance normal to the surface. 


and 


/= [i -«*(-&)] , = ^ 

C^i = 0.081, <T k = 1.0, C e i = 1-44, Cei = 1-83 


O e = 


III RESULTS AND DISCUSSION 

An isotropic eddy-viscosity two-equation model and a nonlinear algebraic stress 
model are evaluated based on their performance in predicting turbulent flow over 
a 13 percent-thick general-aviation airfoil (GA(W)-2). The isotropic eddy-viscosity 
K — e model of Speziale, Abid, and Anderson (SAA) 6 is one of the models; the 
other is the explicit algebraic stress model (EASM) presented in the last section. The 
incompressible = 0.15) experimental data of McGhee et al} is used to assess 
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the models in predicting the aerodynamic characteristics of the airfoil over a range of 
chord Reynolds numbers. 

In the three cases studied (Re = 2.1 x 10 6 , 4.3 x 10®, and 6.3 x 10 6 ), the exper- 
imental transition point is fixed at 7.5 percent chord. In the computations, such a 
forced transition is difficult to implement in an unambiguous manner; however, in the 
computations, because transition occurred near the experimental transition point, no 
additional forcing was necessary. 

Pressure variations over the airfoil at Re = 4.3 x 10 6 at angles of attack of 0° 
and 18° are shown in Figure 1. At a = 0°, both models predict essentially the same 
pressure distribution in the aft portion of the airfoil. Near the leading edge on the 
upper surface, the acceleration of the flow predicted by the K — e model is stronger 
than that predicted by the EASM; and on the lower surface a weaker acceleration 
is predicted by the K — e model. In both cases, the EASM is more accurate in 
predicting the experimental results. At the higher angle of attack, the predictions 
from both models are similar and replicate the experimental results well. 

An important aspect of this study is to assess the ability of the models to correctly 
reproduce the lift and drag characteristics of the airfoils at different angles of attack 
and at different Reynolds numbers. Figure 2 shows the computed lift and drag 
coefficients for various angles of attack (a = 0°, 8°, 14°, 16°, and 18°) at Re = 
2.1 x 10 6 . At this chord Reynolds number, both models correctly predict the stall 
location; however, the EASM performs better than the K — e model in predicting 
the magnitude of Cj; the largest departure from the experimental values occurs at 
the higher angles of attack. This result is in contrast to the previous computational 
studies 2-4 on the ONERA-A airfoil at the same Reynolds number. In these earlier 
studies, the K — e model was unable to accurately predict the correct stall angle, 
and, in fact, did not predict stall at all over the range of angle of attack examined. 
This implies that the details of the airfoil shape are important factors in assessing 
whether a given turbulence model can predict a key aerodynamic characteristic of the 
airfoil. 

The variation of C* with C\ is shown in Figure 2(b). Note that the drag in the 
experiment was determined from a wake- rake (1.17 chords high) survey located 1 
chord length downstream of the airfoil trailing edge. In the computation, the drag 
was deduced from the pressure and viscous force distribution along the airfoil. Nev- 
ertheless, at a = 0°, and to a lesser extent a = 8°, the computational results compare 
reasonably well with experimental results. However, at higher angles of attack, the 
computational results are quantitatively different than those of the experiment. It 
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is not possible to verify whether the wake-rake in the experiment was sufficiently far 
downstream to completely capture the entire wake region at the higher angles of at- 
tack. An analysis of the computed flow field at 1 chord downstream from the trailing 
edge indicated that the rake would not have been sufficient to span the entire wake re- 
gion. This would have led to an underprediction of the drag. Because the experiment 
was performed in a wind tunnel which was approximately 3.75 chords high, while 
the computations were performed in a free-field with computational boundaries at 
15 chords, the possibility axises that in the experiment the rake survey was adequate 
because the wall boundaries would require the flow to turn parallel to the walls and 
thus modify the wake. Unfortunately, it is not possible to conclusively isolate the 
cause of this drag prediction discrepancy at the higher angles of attack within the 
scope of the present study. 

Figure 3 shows the predicted Ci and Cd for Re = 4.3 x 10 6 . At the lower angles of 
attack (a = 0° and 8°), the EASM and to a slightly lesser extent the K — e model, 
does a good job in predicting C\. At higher angles of attack ( a — 16°, 18°, 20°, and 
22°), both models do a poor job in predicting C\. At a = 16° and 18°, both models 
yield the same Ci, but at a = 20° and 22°, the EASM prediction for Ci is less than 
that for the K — e model and is closer to the experimental results. Overall, neither 
model is capable of accurately predicting the C\ near stall. 

As at the lower Reynolds number, Cd is predicted reasonably well at a = 0° and 
8°; however, at the higher angles of attack, the predictions differ significantly. The 
previous discussion on the Cd predictions at the lower Reynolds number is, of course, 
still applicable and would significantly contribute to any discrepancy. 

At the highest Reynolds number studied (Re = 6.3 x 10 6 ), the trends for both Ci 
and Cd are the same as for Re = 4.3 x 10 6 . Figure 4(a) shows that at a = 0° and 8°, 
both models correctly predict the Cf, at the larger angles of attack (a = 15°, 17°, 19°, 
and 21°), neither model accurately predicts the Ci levels. The EASM yields results 
that are closer to the experimental levels than those of the K—e model; however, these 
predictions are still not acceptable. In addition, the stall angle predicted with the 
EASM is shifted slightly from both the K—e model predictions and the experiment al 
value. 

Consistent with results at the other two Reynolds numbers, the comparisons of 
the Cd predictions shown in Figure 4(b) are disappointing. At the higher angles of 
attack, neither model agrees with the experimental Cd results. 

These results indicate that the turbulence models used in this study are incapable 
of predicting the correct airfoil aerodynamic characteristics at higher angles of attack 
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over the Reynolds-number range studied. A clearer picture of this Reynolds-number 
scaling effect is shown in Figure 5 in which the observed experimental behavior can 
be seen to be linear over the range of Reynolds-numbers studied. The K — e model 
displays neither the correct slope nor linear behavior in this region, which suggests 
that a key physical element is not being correctly incorporated into either the calcula- 
tion or the turbulence model. The EASM also predicts the incorrect slope; however, 
the variation across the Reynolds-number range is much closer to linear than for 
the K — e model. Nevertheless, in both computations the qualitative features axe 
predicted inaccurately. 

An explanation for the C/m** predictions can be obtained from an analysis of the 
separation zone at the trailing edge of the airfoil. Unfortunately, the experimental 
study did not include any measurements or visualization of the flow-field near the 
airfoil trailing edge. Figure 6 shows the separation zone predicted by each model along 
the trailing edge of the upper surface near the stall angle of the EASM . At Re = 2.1 x 
10®, the isotropic model does not predict a separation zone along the surface which is 
consistent with the results of Lien and Leschziner 4 ; however, the EASM is found to 
predict separation along the airfoil. When compared to Figure 2, it is seen that the 
EASM performs reasonably well in reproducing both the qualitative and quantitative 
features of the experiment up to the stall angle. As angle of attack increases, both 
models overpredict the Ci value; however, with the EASM predictions, are in better 
agreement with the experimental results. The lack of a predicted separation zone 
along the airfoil by the K — e model suggests that the earlier studies 2-4 at this 
Reynolds number may have completely missed the correct flow behavior along the 
upper surface. 

Figure 6 shows separation zones predicted along the surface by both models at 
the intermediate Reynolds number Re = 4.3 x 10®. Figure 3 shows in this case that 
the K — e model predicts the stall angle more accurately; the EASM overpredicts the 
separation zone at this Reynolds number (Fig. 6). This correlation of separation-zone 
size and lift coefficient suggests that at higher angles of attack both models grossly 
overpredict the size of the separation zone near the trailing edge. At this point in the 
analysis, the results are ambiguous. At Re = 2.1 x 10®, the K — e model does not 
predict separation along the airfoil and, therefore, poorly predicts the C/ distribution 
near stall. However, at Re = 4.3 x 10® the EASM overpredicts the separation-zone 
size and, therefore, underpredicts the C\ distribution near stall. 

The results for Re — 6.3 x 10® are also shown in Figure 6. In this case, the 
EASM more closely approximates the experimental Cj distribution near stall (Fig. 
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4), although the stall angle is underpredicted by both models. Thus, the size of the 
separation zone is underpredicted by the K — e model (Fig. 6). Because both models 
overpredict the C\ distribution at angles of attack greater than the stall angle , it 
can be assumed that at these higher angles of attack the separation-zone size would 
be underpredicted by both models. The reduction in separation-zone size near the 
trailing edge, between the Re = 4.3 x 10 6 and 6.3 x 10® cases, correlates with the 
decrease in shown in Figure 5. 

The results presented in this study have shown that neither the isotropic eddy- 
viscosity K — e model nor the EASM accurately predict the lift and drag coefficients 
over a range of Reynolds numbers and angles of attack. However, areas in which the 
models do perform well relative to one another and in comparison with experiment 
have been identified. Of course, other two-equation K — e models and other ASM’s 
are available for testing; however, while the quantitative results may change slightly, 
the overall qualitative performance characteristics are not expected to change signif- 
icantly. The question then remains to identify which aspect of flow-field physics is 
not being properly modeled. Unfortunately, with the limited amount of mean field 
experimental data available this becomes problematic. 

Although the turbulence models themselves may be inadequate, other dynamic 
features of the flow also may not have been properly accounted for and may need to 
be investigated. One such feature of the experimental flow field, which was not incor- 
porated into the present study, is the unsteady shedding from the trailing edge that 
occurred at and above the intermediate chord Reynolds number (R. J. McGhee, pri- 
vate communication). The present numerical procedure does not solve the governing 
equations in a time-accurate manner and as such cannot give a true representation 
of the time-dependent behavior near the airfoil trailing edge. This unsteadiness may 
affect the separation-zone size, which would result in a different average picture of 
the airfoil characteristics. Some evidence that an important dynamic feature of the 
flow field has been omitted from the calculations can be seen in Figure 5. The lin- 
ear behavior of C/ mar versus Re displayed between Re = 2.1 x 10 6 and 4.3 x 10 6 
(and throughout the Reynolds-number range for the experimental data) is lost at the 
higher Reynolds number, which suggests that some aspect of the physics may have 
been omitted. 

IV CONCLUDING REMARKS 

This purpose of this study was to begin a systematic evaluation of the predictive 
capabilities of two distinct types of turbulence closure models in predicting the aero- 
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dynamic characteristics of airfoils at different Reynolds numbers. One of the models 
was an isotropic eddy- viscosity K — e model; the other was an explicit algebraic stress 
model (EASM). These two models were computationally equivalent in both time and 
memory requirements but were based on significantly different physics. As such, the 
study was not intended to focus on a particular turbulence model within each type 
but on the predictive capabilities of a particular type of model. The selection of 
a different model within each type may have yielded slightly different quantitative 
results, but the overall conclusions would not have differed significantly. 

Neither type of model accurately computed the lift and drag characteristics over 
the Reynolds-number range studied. The EASM more closely predicted the C; dis- 
tributions over a wider range of angle of attack for a particular Reynolds number but 
was unable to accurately predict the distributions near stall. A similar conclusion 
was reached for prediction of the drag coefficient over the Reynolds-number range 
studied, although other factors related to the differences between the computational 
and experimental problems may explain some of these. 

The separation zone at the airfoil trailing edge near stall was also examined at 
the three Reynolds numbers. The predicted separation-zone sizes varied significantly 
among the models; these results were correlated with the lift distributions at the 
three Reynolds numbers. Unfortunately, no unique feature could be identified that 
would indicate the reason for the poor predictive performance of the models at higher 
angles of attack. However, one possible key feature of the experimental study, that 
has not been accurately represented in this study is the unsteady shedding from the 
airfoil trailing edge. Shedding would affect the separation-zone size predicted by the 
different models and could lead to different average values for Cj and C&. The next 
phase of this study will address this issue by performing these same calculations with 
a time-accurate numerical solver. 
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Figure 1. Pressure variation around airfoil at a = 0° and 18° with Re = 4.3 x 10 6 . 
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Figure 2. Lift and drag airfoil characteristics at Re = 2.1 x 10®. (a) Variation of lift 
coefficient with angle of attack, (b) Variation of drag and lift coefficients. 
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Figure 4. Lift and drag airfoil characteristics at Re = 6.3 x 10®. (a) Variation of lift 
coefficient with angle of attack, (b) Variation of drag and lift coefficients. 
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Figure 5. Variation of maximum lift coefficient with Reynolds number. 
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